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Abstract 



Recently we presented a simple method for determining the correlated un- 
certainties of the light element abundances expected from big bang nucle- 
osynthesis, which avoids the need for lengthy Monte Carlo simulations. We 
now extend this approach to consider departures from the Standard Model, 
in particular to constrain any new light degrees of freedom present in the 
thermal plasma during nucleosynthesis. Since the observational situation re- 
garding the inferred primordial abundances has not yet stabilized, we present 
illustrative bounds on the equivalent number of neutrino species N v for var- 
ious combinations of individual abundance determinations. Our 95% C.L. 
bounds on N u range between 2 and 4, and can easily be reevaluated using the 
technique provided when the abundances are known more accurately. 
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I. INTRODUCTION 



The Standard Model (SM) contains only N v = 3 weakly interacting massless neutrinos. 
However the recent experimental evidence for neutrino oscillations [0J may require it to be 
extended to include new superweakly interacting massless (or very light) particles such as 
singlet neutrinos or Majorons. These do not couple to the Z° vector boson and are therefore 
not constrained by the precision studies of Z° decays which establish the number of <S77(2)l 
doublet neutrino species to be || 

N v = 2.993 ±0.011. (1) 

However, as was emphasized some time ago ||, such particles would boost the relativistic 
energy density, hence the expansion rate, during big bang nucleosynthesis (BBN), thus 
increasing the yield of 4 He. This argument was quantified for new types of neutrinos and 
new superweakly interacting particles || in terms of a bound on the equivalent number of 
massless neutrinos present during nucleosynthesis: 

iV„ = 3 + /B, F Ef (^T) 4 , (2) 



- V 



where <?j is the number of (interacting) helicity states, /b = 8/7 (bosons) and /p = 1 
(fermions), and the ratio Ti/T v depends on the thermal history of the particle under consid- 
eration ||. For example, Ti/T v < 0.465 for a particle which decouples above the electroweak 
scale such as a singlet Majoron or a sterile neutrino. However the situation may be more 
complicated, e.g. if the sterile neutrino has large mixing with a left-handed doublet species, 
it can be brought into equilibrium through (matter-enhanced) oscillations in the early uni- 
verse, making TijT v ~ 1 ||. Moreover such oscillations can generate an asymmetry between 
u e and z/ e , thus directly affecting neutron-proton interconversions and the resultant yield of 
4 He 0. This can be quantified in terms of the effective value of N v parametrizing the expan- 
sion rate during BBN, which may well be below 3! Similarly, non-trivial changes in N u can 
be induced by the decays or annihilations || of massive neutrinos (into e.g. Majorons), 
so it is clear that it is a sensitive probe of new physics. 

The precise bound on N v from nucleosynthesis depends on the adopted primordial ele- 
mental abundances as well as uncertainties in the predicted values. Although the theoretical 
calculation of the primordial 4 He abundance is now believed to be accurate to within ±0.4% 
HI , its observationally inferred value as reported by different groups differs by as 



much as ~ 4%. Furthermore, a bound on N v can only be derived if the nucleon-to-photon 
ratio 1] = n-^/n^ (or its lower bound) is known, since the effect of a faster expansion rate 
can be compensated for by the effect of a smaller nucleon density. This involves comparison 
of the expected and observed abundances of other elements such as D, 3 He and 7 Li which 
are much more poorly determined, both observationally and theoretically. The most crucial 
element in this context is deuterium which is supposedly always destroyed and never cre- 
ated by astrophysical processes following the big bang |TB[ . Until relatively recently TE 



its primordial abundance could not be directly measured and only an indirect upper limit 
could be derived based on models of galactic chemical evolution. As reviewed in ref . , the 
implied lower bound to r\ was then used to set increasingly stringent upper bounds on N v 
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ranging from 4 downwards |17j , culminating in one below 3 which precipitated the so-called 
for standard BBN 



crisis 



18 



However as cautioned before UI9 



and was interpreted as requiring new physics. 

there are large systematic uncertainties in such con- 



straints on N v which are sensitive to our limited understanding of galactic chemical evolu- 
tion. Moreover it was emphasized [[2(J that the procedure used earlier |I7j to bound N v was 
statistically inconsistent since, e.g., correlations between the different elemental abundances 
were not taken into account. A Monte Carlo (MC) method was developed for estimation 
of the correlated uncertainties in the abundances of the synthesized elements [21, 22], and 
incorporated into the standard BBN computer code |J23[| , thus permitting reliable determi- 
nation of the bound on N u from estimates of the primordial elemental abundances. Using 
this method, it was shown p4f that the conservative observational limits on the primordial 
abundances of D, 4 He and 7 Li allowed N u < 4.53 (95% C.L.), significantly less restrictive 
than earlier estimates. Similar conclusions followed from studies using maximum likelihood 
(ML) methods PoVf27|j. However the use of the Monte Carlo method is computationally 
expensive and moreover the calculations need to be repeated whenever any of the input 
parameters — either reaction rates or inferred primordial abundances — are updated. 

In a previous paper [28] we presented a simple method for estimation of the BBN abun- 
dance uncertainties and their correlations, based on linear error propagation. To illustrate 
its advantages over the MC+ML method, we used simple \ 2 statistics to obtain the best-fit 
value of the nucleon-to-photon ratio in the standard BBN model with N v = 3 and indi- 
cated the relative importance of different nuclear reactions in determining the synthesized 
abundances. In this work we extend this approach to consider departures from N u = 3. We 
have checked that our results are consistent with those obtained independently [^9] using 
the MC+ML method [29,30] where comparison is possible. 

The essential advantage of our method is that the correlated constraints on N u and r\ 
can be easily reevaluated using just a pocket calculator and the numerical tables provided, 
when the input nuclear reaction cross-sections or inferred abundances are known better. We 
have in fact embedded the calculations in a compact Fortran code, which is available upon 
request from the authors, or from a website |H[ . Thus observers will be able to readily assess 
the impact of new elemental abundance determinations on an important probe of physics 
beyond the standard model. 



II. THE METHOD 



In this section we recapitulate the basics of our method |28[ and outline its extension to 
the case N v ^ 3. 



A. Basic Ingredients 



The method has both experimental and theoretical ingredients. The experimental in- 
gredients are: (a) the inferred values of the primordial abundances, Yi ± af, and (b) the 
nuclear reaction rates, Rk ± Ai?^. We normalize [28 
values (Rk = 1), namely, to the values compiled in Ref. []2"2 
rate, which is updated to its current value fl2j. 



i 3 

all the rates to a "default" set of 
except for the neutron decay 
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The theoretical ingredients are: (a) the calculated abundances Yf, and (b) the logarithmic 
derivatives Xik = <9 In 1^/(9 In R^. Such functions have to be calculated for generic values of 
N u and of x = log 10 (r] 10 ), where i] 10 = r]/10~ w . Note that the fraction of the critical density 
in nucleons is given by Q^h 2 — ^io/273, where h ~ 0.7±0.1 is the present Hubble parameter 
in units of 100 km s _1 Mpc -1 , and the present temperature of the relic radiation background 
is T = 2.728 ±0.002 K @. 

The logarithmic derivatives Xik can be used to to propagate possible changes or updates 
of the input reaction rates (Rk — > Rk + SRk) to the theoretical abundances (Yj — > Yj + 
Yi\ik5 Rk I Rk) • Moreover, they enter in the calculation of the theoretical error matrix for 
the abundances, cr^ = Yj Yj J2k \k^jk(^ Rk/ Rk) 2 - This matrix, summed to the experimental 
error matrix afj = SijaiWj and then inverted |28f , defines the covariance matrix of a simple x 2 
statistical estimator. Contours of equal x 2 can then be used to set bounds on the parameters 
(x, N v ) at selected confidence levels. 

In Ref. [^8] we gave polynomial fits for the functions Yi(x, N v ) and Xik(x, N u ) for x G [0, 1] 
and N v = 3. The extension of our method to the case N v ^ 3 (say, 1 < N v < 5) is, in 
principle, straightforward, since it simply requires recalculation of the functions Y^ and 
Xik at the chosen value of iV,,. However, it would not be practical to present, or to use, 
extensive tables of polynomial coefficients for many different values of N u . Therefore, we 
have devised some formulae which, to good accuracy, relate the calculations for arbitrary 
values of N u to the standard case N v = 3, thus reducing the numerical task dramatically. 
Such approximations are discussed in the next subsection. 



B. Useful Approximations 

As is known from previous work ||32|| , the synthesized elemental abundances D/H, 3 He/H, 
and 7 Li/H (i.e., Y 2 , Y 3 , and Y? in our notation) are given to a good approximation by the 
quasi-fixed points of the corresponding rate equations, which formally read 

dY 

^«^y x y X H Tl (3) 

where the sum runs over the relevant source (+) and sink (— ) terms, and (ctv)t is 
the thermally-averaged reaction cross section. Since the temperature of the universe 
evolves as dT/dt oc — T 3 ^/^", with the number of relativistic degrees of freedom, g+ = 
2 + (7/4)(4/ll) 4/ ' 3 A ;/ (following e + e~ annihilation), the above equation can be rewritten as 

^^-^r 2 T- 3 T,YxYx(av) T , (4) 
ai 9* +,- 

which shows that Y 2 , Y 3 , and Y 7 depend on rj and N u essentially through the combina- 

1 /2 

tion rj/gj . Thus the calculated abundances Y 2 , I3, and Y 7 (as well as their logarithmic 
derivatives Xik) should be approximately constant for 



1 

log 77 - - log^ = const , (5) 



as we have verified numerically. 
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Equation (|), linearized, suggests that the values of Yi and of \ k for N u = 3 + AN U can 
be related to the case N v = 3 through an appropriate shift in x: 

Yi{x, 3 + AN V ) ~Yi{x + Ci AN„, 3) , (6) 
X ik {x, 3 + AN„) ~ \ ik (x + c { AN u , 3) , (7) 

where the coefficient c, is estimated to be ~ —0.03 from Eq. (^) (at least for small AN U ). In 
order to obtain a satisfactory accuracy in the whole range (x, N v ) G [0, 1] x [1, 5], we allow 
upto a second-order variation in AN U , and for a rescaling factor of the Y^s: 

Yi(x, 3 + AN U ) = (1 + ai AN v + biAN u 2 ) Y { {x + Ci AN u + diANj 1 , 3) , (8) 
X ik (x, 3 + AN U ) = X ik {x + qAN u + diAN u 2 , 3) . (9) 

We have checked that the above formulae (with coefficients determined through a numerical 
best-fit) link the cases N v ^ 3 to the standard case N u = 3 with very good accuracy. 

As regards the 4 He abundance (Y4 in our notation), a semi-analytical approximation also 
suggests a relation between x and N v similar to Eq. ((^), although with different coefficients 
PB| . Indeed, functional relations of the kind work well also in this case. However, in 
order to achieve higher accuracy and, in particular, to match the result of the recent precision 
calculation of Y4 which includes all finite temperature and finite density corrections [ID|, we 
also allow for a rescaling factor for the A^'s. 

We wish to emphasize that the validity of our prescription for the evaluation of the 
BBN uncertainties and for the x 2 statistical analysis does not depend on the approximations 
discussed above. The semi-empirical relations (|8||]) are only used to enable us to provide the 



interested reader with a simple and compact numerical code |31| . This allows easy extraction 
of joint fits to x and N v for a given set of elemental abundances, without having to run the 
full BBN code, and with no significant loss in accuracy. 



III. PRIMORDIAL LIGHT ELEMENT ABUNDANCES 

The abundances of the light elements synthesized in the big bang have been subsequently 
modified through chemical evolution of the astrophysical environments where they are mea- 
sured The observational strategy then is to identify sites which have undergone as little 



chemical processing as possible and rely on empirical methods to infer the primordial abun- 
dance. For example, measurements of deuterium (D) can now be made in quasar absorption 
line systems (QAS) at high red shift; if there is a "ceiling" to the abundance in different 
QAS then it can be assumed to be the primordial value. The helium ( 4 He) abundance is 
measured in H II regions in blue compact galaxies (BCGs) which have undergone very lit- 
tle star formation; its primordial value is inferred either by using the associated nitrogen 
or oxygen abundance to track the stellar production of helium, or by simply observing the 
most metal-poor objects J35j. (We do not consider 3 He which can undergo both creation and 
destruction in stars [|33J and is thus unreliable for use as a cosmological probe.) Closer to 
home, the observed uniform abundance of lithium ( 7 Li) in the hottest and most metal-poor 
Pop II stars in our Galaxy is believed to reflect its primordial value . 



However as observational methods have become more sophisticated, the situation has 
become more, instead of less, uncertain. Large discrepancies, of a systematic nature, have 
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emerged between different observers who report, e.g., relatively 'high' or 'low 



15,39,40|| values of deuterium in different QAS, and 'low' [llj4l| or 'high' [12,42] values 



of helium in BCG, using different data reduction methods. It has been argued that the 
Pop II lithium abundance |43| [45| may in fact have been significantly depleted down from 

We do not wish to 



its primordial value fl4"6"||4"Tf , with observers arguing to the contrary |4£ 
take sides in this matter and instead consider several combinations of observational deter- 
minations, which cover a wide range of possibilities, in order to demonstrate our method 
and obtain illustrative best-fits for r] and N u . The reader is invited to use the programme 
we have provided JTIJ to analyse other possible combinations of observational data. 



A. Data Sets 



The data sets we consider are tabulated in Table |. Below we comment in detail on our 
choices. 



Data Set A: This is taken from Ref. |29| who performed the first detailed MC+ML 
analysis to determine 77 and N v and is chosen essentially for comparison with our 
method, as in our previous paper pf| . 

Their adopted value of the primordial deuterium abundance, Y 2 = 1.9 ±0.4 x 10" 4 , was 
based on early observations of a QAS at redshift z = 3.32 towards Q0014+813 which 
suggested a relatively 'high' value ||14j| , and was consistent with limits set in other 
QAS, but in conflict with the much lower abundance found in a QAS at z = 3.572 
towards Q1937-1009 |15| . More recently, observations of a QAS at z = 0.701 towards 



Q1718+4807 have also yielded a high abundance PT| , |3"8"|j as we discuss later. 

The primordial helium abundance was taken to be I4 = 0.234±0.002±0.005 from linear 
regression to zero metallicity in a set of 62 BCGs [£|l|], based largely on observations 



which gave a relatively 'low' value |11 



Finally the primordial lithium abundance Y 7 = 1.6 ± 0.36 x 10 10 was taken from the 
Pop II observations of Ref. |[44||, assuming no depletion. 



Data Set B: This corresponds to the alternative combination of 'low' deuterium and 



'high' helium, as considered in our previous work ||28|| , with some small changes. 

The primordial deuterium abundance, Y 2 = 3.4±0.3x 10~ 5 , adopted here is the average 
of the 'low' values found in two well-observed QAS, at z = 3.572 towards Q1937-1009 
0, and at z = 2.504 towards Q1009+2956 0. 



The primordial helium abundance, Y4 = 0.245 ± 0.004, is taken to be the average of 
the values found in the two most metal-poor BCGs, I Zw 18 and SBS 0335-052, from a 
new analysis which uses the helium lines themselves to self-consistently determine the 
physical conditions in the H II region, and specifically excludes those regions which are 



believed to be affected by underlying stellar absorption p2| . For example these authors 
demonstrate that there is strong underlying stellar absorption in the NW component 
of I Zw 18, which has been included in earlier analyses [p| . 
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The primordial lithium abundance Yy = 1.73 ± 0.21 x 10 10 is from Ref. ||45|| , again 



assuming no depletion. (Note that the uncertainty was incorrectly reported as ±0.12 x 
10~ 10 in Ref. |j36| , as used in our previous work p8| .) 



Data Set C: It has been suggested [|49J that the discordance between the 'high' and 
'low' values of the deuterium abundance reported in QAS may be considerably re- 
duced if the analysis of the H+D profiles accounts for the correlated velocity field of 
bulk motion, i.e. mesoturbulence, rather than being based on multi-component mi- 
croturbulent models. It is then found |49| that a value of Y 2 = (3.5 — 5.2) x 10~ 5 is 
compatible simultaneously (at 95% C.L.) with observations of the QAS at z — 0.701 
towards Q1718+4807 (in which a 'high' value was reported fl37,38|), and observations 



of the QAS at z = 3.572 towards Q1937-1009 and at z = 2.504 towards Q1009+2956 
(in which a 'low' value was reported |39| , [40| ). We adopt this value, along with the 
same helium abundance as in Set B. 

It has also been argued that the lithium abundance observed in Pop II stars has 
been depleted down from a primordial value of Y? = 3.9 ± 0.85 x 10~ 10 [50|, the 



lower end of the range being set by the presence of highly overdepleted halo stars 
and consistency with the 7 Li abundance in the Sun and in open clusters, while the 
upper end of the range is set by the observed dispersion of the Pop II abundance 
"plateau" and the 6 Li/ 7 Li depletion ratio. We adopt this value, noting that a somewhat 
smaller depletion is suggested by other workers [[171] who find a primordial abundance 
of Y 7 = 2.3 ±0.5 x 10- 10 . 

Data Set D: Recently, a 'high' value of the deuterium abundance, Y% = 3.3±1.2 x 10~ 4 , 



has been reported from observations of a QAS at z = 0.701 towards Q1718±4807 |38 
in confirmation of an earlier claim |37j. We adopt this value along with the same 
helium abundance as in set A. 

For the lithium abundance, we adopt the same value jy| as in Set B but increase the 



systematic error by 0.02 dex to allow for the uncertainty in the oscillator strengths of 
the lithium lines |5l|]. 



B. Qualitative Implications on N v and r\ 

Different choices for the input data sets (A-D) lead to different implications for rj and 
N u , that can be qualitatively understood through Figs. |l]-4l[ 

Figure |1] shows the BBN primordial abundances Yi (solid lines) and their ±2a bands 
(dashed lines), as functions of x = log(^/10 -10 ), for N v = 2, 3, and 4. The grey areas 
represent the ±2<r bands allowed by the data set A (see Table |). There is global consistency 
between theory and data for x ~ 0.2 — 0.4 and N v = 3, while for N v = 2 (N u = 4) the Y2 
data prefer values of x lower (higher) than the Y4 data. Therefore, we expect that a global 
fit will favor values of (x, N v ) close to (0.3, 3). 

Figure |] is analogous, but for the data set B. In this case, there is still consistency 
between theory and data at N v = 3, although for values of x higher than for the data set 
A. For N v = 2 (N u = 4) the combination of Y 2 and Y? data prefer values of x lower (higher) 
than Y4. The best fit is thus expected to be around (x, iVj,) ~ (0.7, 3). 
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Similarly, Figure |] shows the abundances for the data set C. The situation is similar to 
data set B (Fig. §), but one can envisage a best fit at a slightly lower value of x, due to the 
higher value of Y 2 , partly opposed by the increase in Y 7 . 

Finally, Fig. ^ refers to the data set D. In this case, data and theory are not consistent 
for N v = 2, since Y 2 and Y4 pull x in different directions, and no "compromise" is possible 
since intermediate values of x are disfavored by Y 7 . However, for N v = 3 there is relatively 
good agreement between data and theory at low x. Therefore, we expect a best fit around 
(x,N v )~ (0.2,3). 

The qualitative indications discussed here are confirmed by a more accurate analysis, 
whose results are reported in the next section. 



IV. DETERMINING N v 

We now present the results of fits to the data sets A-D in the (x, N v ) variables, using our 
method to estimate the correlated theoretical uncertainties, and adopting \ 2 statistics to 
include both theoretical and experimental errors. We have used the theoretical Yj's obtained 



by the standard (updated) BBN evolution code p3| , and checked that using the polynomial 
fits given in Sec. II B induce negligible changes which would not be noticeable on the plots. 

In the analysis, we optionally take into account a further constraint on rj (independent on 
N v ) coming from a recent analysis of the Lya- "forest" absorption lines in quasar spectra. The 
observed mean opacity of the lines requires some minimum amount of neutral hydrogen in 
the high redshift intergalactic medium, given a lower bound to the flux of ionizing radiation. 
Taking the latter from a conservative estimate of the integrated UV background due to 
quasars, Ref. finds the constraint rj > 3.4 x 10~ 10 . This bound is not well-defined 
statistically but, for the sake of illustration, we have parametrized it through a penalty 
function quadratic in rj: 

xu,h 2.7 x (Mi^y , (10) 



V 



to be eventually added to the X 2 (.ViN u ) derived from our fit to the elemental abundances. 
The above function excludes values of rj smaller than 3.4 x 10~ 10 at 90% CL. (for one degree 
of freedom, rj). 

Figure [| shows the results of joints fits to x — log(r/io) an d N v using the abundances of 
data set A. The abundances Y 2 , I4, and Y 7 are used separately (upper panels), in combina- 
tions of two (middle panels), and all together, without and with the Lya-forest constraint 
on rj (lower panels). In this way the relative weight of each piece of data in the global fit can 
be understood at glance. The three CL. curves (solid, thick solid, and dashed) are defined 
by X 2 ~ Xmin = 2-3, 6.2, and 11.8, respectively, corresponding to 68.3%, 95.4%, and 99.7% 
CL. for two degrees of freedom (77 and N u ), i.e., to the probability intervals often designated 
as 1, 2, and 3 standard deviation limits. The \ 2 is minimized for each combination of Y it 
but the actual value of Xmin (and the best-fit point) is shown only for the relevant global 
combination Y 2 + I4 + Y 7 (+Lya). 

The results shown in Fig. [|for the combinations Y4+Y7 and Y 2 -\-Y±-\-Y 7 are consistent with 
those obtained in ref. |29j by using the same input data but a completely different analysis 
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method (namely Monte Carlo simulation plus maximum likelihood). The consistency is 
reassuring and confirms the validity of our method. For this data set, the helium and 
deuterium abundances dominate the fit, as it can be seen by comparing the combinations 
Y 2 + Y4 and Y 2 + 14 + Y 7 . The preferred values of x are relatively low, and the preferred 
values of N u range between 2 and 4. Although the fit is excellent, the low value of x is in 
conflict with the Lya-forest constraint on rj, as indicated by the increase of x min from 0.02 
to 8.89. 

Figure |6] is analogous, but for the data set B which favors high values of x because of 
the 'low' deuterium abundance. The combination of Y 2 + Y 7 isolates, at high x, a narrow 
strip which depends mildly on N u . The inclusion of Y4 selects the central part of such strip, 
corresponding to N v between 2 and 4. As in Fig. |5|, the combination Y±-\-Y 7 does not appear 
to be very constraining. The overall fit to Y 2 + Y4 + Y 7 is acceptable but not particularly 
good, mainly because Y 2 and Y 7 are only marginally compatible at high x. On the other 
hand, the Y 2 + Y4 + Y 7 bounds are quite consistent with the Lya-forest constraint. 

In data set C, the deuterium abundance has increased further. Moreover the lithium 
abundance is no longer at the minimum of the theoretical curve as before, so strongly 
disfavors "intermediate" values of x. The overall effect, as shown in Figure [7|, is that Xmm 
decreases a bit with respect to Set B, and the best fit value of x moves slightly lower. The 
allowed value of N u ranges between 2 and 4. Note that had we retained the same Y 7 as in 
Set B, then Xmm would have dropped to 0.91 (2.55) for the combination Y 2 + Y A + Y 7 (+ 
Lya-forest constraint). 

Finally, Fig. |8| refers to data set D which, like Data Set A, has the 'high' deuterium 
abundance but with larger uncertainties. So although a low value of x is still picked out, a 
high x region is still possible at the 2a level (in the Y 2 + Y 4 + Y 7 panel) and is even favored 
when the Lya-forest bound is included (although with an unacceptably high Xmm)- Note 
that had the lithium abundance been taken to be the same as in data set C (i.e. allowing 
for depletion), the Xmm would have been 0.07 (7.42) for the combination Y 2 + Y4 + Y 7 (+ 
Lya-forest constraint). 

Of course one can also consider orthogonal combinations to those above, e.g. 'high' deu- 
terium and 'high' helium, or 'low' deuterium and 'low' helium ||30|| . The latter combination 
implies N u ~ 2, thus creating the so-called "crisis" for standard nucleosynthesis Con- 



versely, the former combination suggests N v ~ 4, which would also constitute evidence for 
new physics. Allowing for depletion of the primordial lithium abundance to its Pop II value, 
relaxes the upper bound on N v further, as noted earlier \2~A\ . 



V. CONCLUSIONS 

The results discussed above demonstrate that the present observational data on the 
primordial elemental abundances are not as yet sufficiently stable to derive firm bounds on 
r\ and N v . Different and arguably equally acceptable choices for the input data sets lead 
to very different predictions for rj, and to relatively loose constraints on N v in the range 
2 to 4 at the 95% C.L. Thus it may be premature to quote restrictive bounds based on 
some particular combination of the observations, until the discrepancies between different 
estimates are satisfactorily resolved. Our method of analysis provides the reader with an 
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easy-to-use technique |31] to recalculate the best-fit values as the observational situation 
evolves further. 

However one might ask what would happen if these discrepancies remain? We have 
already noted the importance of an independent constraint on 77 (from the Lya-forest) in 
discriminating between different options. However, given the many assumptions which go 
into the argument |5^| , this constraint is rather uncertain at present. Fortunately it should 
be possible in the near future to independently determine rj to within ~ 5% through mea- 
surements of the angular anisotropy of the cosmic microwave background (CMB) on small 
angular scales |53], in particular with data from the all-sky surveyors MAP and PLANCK 



54|. Such observations will also provide a precision measure of the relativistic particle 
content of the primordial plasma. Hopefully the primordial abundance of 4 He would have 
stabilized by then, thus providing, in conjunction with the above measurements, a reliable 
probe of a wide variety of new physics which can affect nucleosynthesis. 
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TABLES 



TABLE I. Experimental data sets considered in this paper for the elemental abundances Y\. 





A 


B 


C 


D 


Y 2 x 10 5 


19 ±4 


3.4 ±0.3 


4.35 ± 0.43 


33 ± 12 


T 4 


0.234 ± 0.0054 


0.245 ± 0.004 


0.245 ± 0.004 


0.234 ± 0.0054 


Y 7 x 10 10 


1.6 ±0.36 


1.73 ±0.21 


3.9 ±0.85 


1.73 ±0.29 
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FIGURES 



Data set ! 

Y 2 = (1 .9 ± 0.4) x 10" 4 
Y 4 = 0.234 ± 0.0054 
Y 7 = (1.6 ± 0.36) x 10~ 10 

N„ = 2 



0.26 



I 0.24 - 



- 0.22 - 



N„ = 3 





Allowed bands '. 

^±±^±±^ THEORY (± 2a) 
DATA (±2ct) 

N„ = 4 



0.26 



- 0.24 



- 0.22 - 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 





0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 



i°g(7?io) 



i°g(?7io) 



FIG. 1. Primordial abundances Y4 ( 4 He mass fraction) Y2 (D/H) and Yy ( 7 Li/H), for N v = 2, 
3, and 4. Solid and dashed curves represent the theoretical central values and the ±2<r bands, 
respectively. The grey areas represent the ±2<r experimental bands for the data set A in Table |. 
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Data set 



Y 2 = (3.4 ± 0.3) x 10" 5 

Y 4 = 0.245 ± 0.004 

Y 7 = (1 .73 ± 0.21) x 10 -10 



Allowed bands '. 

THEORY (± 2a) 
DATA (± 2a) 
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Data set 



Y 2 = (4.35 ± 0.43) x 1 0" 5 

Y 4 = 0.245 ± 0.004 

Y 7 = (3.9 ± 0.85) x 10~ 10 



Allowed bands '. 

THEORY (± 2a) 
DATA (± 2a) 
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Data set 



Y 2 = (3.3 ± 1 .2) x 10" 4 

Y 4 = 0.234 ± 0.0054 

Y 7 = (1 .73 ± 0.29) x 10 -10 



Allowed bands '. 

THEORY (± 2a) 
DATA (± 2a) 
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Data set 

Y 2 = (1.9 ± 0.4) x 10" 4 
Y 4 = 0.234 ± 0.0054 
Y 7 = (1.6 ± 0.36) x 10" 



Allowed regions ! 



68.3% C.L 
95.4% C.L. 
99.7% C.L. 



N, 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 



N, 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



N, 



3 - 



2 - 




0.2 0.4 0.6 



1 0.2 0.4 0.6 0. 

Iog(??io) 



1 0.2 0.4 0.6 0. 



Iog(?7 10 ) 



FIG. 5. Joints fits to x = log(r/io) and N u using the abundances of data set A. The abundances 
Yz, Y4, and Y? are used separately (upper panels), in combinations of two (middle panels), and all 
together, without and with the Lya-forest constraint on 77 (lower panels). 



18 



Data set ! 

Y 2 = (3.4 ± 0.3) x 10" 5 

Y 4 = 0.245 ± 0.004 

Y 7 = (1 .73 ± 0.21) x 10" 10 



Allowed regions ! 



68.3% C.L 
95.4% C.L. 
99.7% C.L. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



N, 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



5 




5 








5 






; y 2 +y 4 +y 7 ; 




- Lya 








; Y 2 +Y 4 +Y 7 +Lya 


4 




4 










4 




N, 


















3 




3 










3 




2 




2 










2 






- X min 3.71 














- X min 4.94 


1 




1 


I, , 








1 





0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

log(77 10 ) log(?7 10 ) log(^ 10 ) 



FIG. 6. As in Fig. ||, but for the data set B. 
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Data set ! 

Y 2 = (4.35 ± 0.43) x 10" 

Y 4 = 0.245 ± 0.004 

Y 7 = (3.90 ± 0.85) x 10" 



Allowed regions ! 



68.3% C.L 
95.4% C.L. 
99.7% C.L. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
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: xL = 2.80 ; 
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I, , 








1 


: y 2 - = 4 1 7 

A mm ** 1 ' 



0.2 0.4 0.6 0. 

Iog(r?io) 



1 0.2 0.4 0.6 0. 

Iog(7? 10 ) 



1 0.2 0.4 0.6 0. 

Iog(>?io) 



FIG. 7. As in Fig. ||, but for the data set C. 
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Data set ! 

Y 2 = (3.3 ± 1.2) x 10" 4 
Y 4 = 0.234 ± 0.0054 
Y 7 = (1 .73 ± 0.29) x 10" 



Allowed regions ! 



68.3% C.L 
95.4% C.L. 
99.7% C.L. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



N, 



3 

2 

1 t 



'Y 2 +Y 4 










0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



N, 



2 - 




l 5 


- Lya 






: 4 










: 3 










: 2 










J 1 


I, , 









Y 2 +Y 4 +Y 7 +Lya 




X min 7.59. 



0.2 0.4 0.6 0. 

Iog(r?io) 



log(7? 10 ) 



log(7?io) 



1 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



FIG. 8. As in Fig. §, but for the data set D. 
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